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摘 要 目前 调节 效应 检验 主要 是 基于 截面 数据 ,本 文 讨论 纵向 (追踪 ) 数 据 的 调节 效应 分 析 。 如 果 自 变量 X 
和 因 变 量 了 有 纵向 数据 ,调节 效应 可 分 为 三 类 : 调节 变量 Z 不 随时 间 变 化 、Z 随时 间 变 化 、 调 节 变 量 从 自 变 
量 或 因 交 量 中 产生 。 评 介 了 基于 多 层 模 型 、 多 层 结构 方程 模型 、 交 又 滞后 模型 和 潜 交 量 增 长 模型 的 纵向 数据 
的 多 种 调节 效应 分 析 方法 。 调 节 效 应 的 分 解 和 潜 调 节 结 构 方程 法 的 使 用 是 纵向 数据 的 调节 效应 分 析 的 两 大 特 
点 。 对 基于 四 类 模型 的 调节 效应 分 析 方 法 进行 综合 比较 后 ， 总 结 出 一 个 纵向 数据 的 调节 效应 分 析 流 程 。 随 后 
用 实际 例子 演示 如 何 进 行 纵 向 数据 的 调节 效应 分 析 ， 并 给 出 相应 的 Mplus 程序 。 随 后 展望 了 纵向 数据 的 调节 
效应 分 析 的 拓展 方向 , 例如 基于 动态 结构 方程 模型 的 密集 追踪 数据 的 调节 效应 分 析 。 

关键 词 ” 纵 向 数据 ,调节 效应 ， 多 层 模 型 ， 多 层 结构 方程 模型 交 又 灌 后 模型 ， 潜 交 量 增长 模型 
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调节 (moderatiom) 是 社会 科学 研究 中 重要 的 增长 模型 (Latent Growth Model, LGM) 的 纵向 数据 
方法 学 概念 。 如 果 自 变量 邓 对 因 变 量 了 的 作用 随 的 调节 效应 分 析 方 法 后 ， 对 基于 4 类 模型 的 调节 
第 三 个 变量 Z 的 变化 而 变化 , 则 称 Z 在 和 和 了 之 效应 分 析 方 法 进行 了 综合 比较 , 并 提出 一 套 纵向 
间 起 调节 作用 ， 此 时 称 Z 为 调节 变量 ( 见 图 1)。 近 数据 的 调节 效应 分 析 流程 ; 然后 用 实际 例子 演示 
30 年 来 ， 调 节 研 究 备 受 关注 ， 分 析 方 法 和 实际 应 如 何 进行 纵向 数据 的 调节 效应 分 析 。 最 后 展望 了 
用 都 得 到 长 足 发 展 。 然 而 目前 调节 效应 分 析 主 要 纵向 数据 的 调节 效应 分 析 的 拓展 方向 。 
涉及 截面 数据 (cross-sectional data), 许多 研究 者 
对 截面 数据 的 调节 效应 分 析 早 已 驾轻就熟 ， 但 是 ， 
当面 对 纵向 (追踪 ) 数 据 (longitudinal data) 时 ， 如 何 
进行 调节 效应 分 析 呢 ? 近年 来 ， 纵向 数据 的 调节 
效应 分 析 已 经 成 为 调节 研究 的 新 热点 ( 方 杰 等 ， 
2018; GA Be, 2017; W BE, XY ZL zs, 2020; 


Ozkok et al., in press; Preacher et al., 2016; Wen et (a) 未 中 心 化 (b) 中 心 化 
ee 图 1 基于 多 层 模型 的 纵向 数据 的 调节 模型 


Multilevel Modeling,，MLM)、 多 层 结构 方程 模型 


(Multilevel Structural Equation Model, MSEM)、 交 1 基于 多 层 模 型 的 纵向 数据 的 调节 效 
又 滞后 模型 (Cross-Lagged Model, CLM) 和 港 变量 应 分 析 


纵向 (追踪 ) 数 据 可 看 成 是 一 个 两 层 结构 ， 重 
收 稿 日 期 : 2021-10-16 复 测量 或 测量 点 1 为 第 1 层 ,， 观测 个 体 j 为 第 2 层 ， 
* 国家 自然 科学 基金 项 目 (32171091)、 国 家 社会 科学 基 可 使 用 多 层 模 型 进行 纵向 数据 分 析 。 如 果 调 节 变 


金 项 目 (17BTJ035) 资 助 。 Jra 
通信 作者 : aL BE wenzl@scnu.edu.cn 量 不 随时 间 变 化 (即时 不 变 ， 只 需要 一 次 测量 )， 


at 


2461 


iv 合作 期 刊 


2462 心理 科学 进展 


ChinaXiv 合 作 期 干 


第 30% 


则 调节 变量 在 层 2 (用 Zz 表示 )， 自 变量 和 因 变 量 
都 在 层 1， 可 使 用 2x(1 一 D) 调 节 模 型 (这 三 个 数字 
依次 代表 调节 变量 、 自 变量 和 因 变 量 的 层级 ， 数 
字 2 表 示 层 级 2, 数字 1 表示 层级 1) 进 行 纵向 数据 
的 调节 效应 分 析 ， 具体 分 析 方 法 详 见 综述 方 杰 等 
(2018), 不 再 袭 述 。 但 如 果 调 节 变 量 是 随时 间 变 化 
(即时 变 ) 的 ， 可 以 像 自 变量 和 因 变 量 那 样 获 取 纵 
向 数据 (用 Zy 表示 )， 此 时 的 调节 效应 分 析 复 杂 多 
了 ， 以 下 将 详 述 相应 的 调节 效应 分 析 方 法 。 
1.1 传统 的 多 层 调 节 效应 分 析 

如 果 所 有 变量 都 随时 间 变 化 , 则 所 有 变量 都 
在 层 1, 可 使 用 1x(1 一 1) 调 节 模 型 进行 纵向 数据 
的 调节 效应 分 析 。 此 时 , 纵向 数据 的 调节 效应 分 
析 ( 图 1(a), Preacher et al., 2016) 可 表示 为 

层 -1 : =poj+ BjXy+ 


By Zu + Bs, XyZy + Ey 0) 

层 -2: Boj =o + Mo; (2a) 
Bip =Not lh; Poy =+ Lb; 

Ps; = 39+ Hs; (2b) 


将 方程 (2a) 和 (2b) 代 入 方程 (1) 得 : 
Yi = o0 + 710X y + ¥20Zy + V30X GZy + (3) 
My jXq + Ma jbg + taj XqZy + Ey + Maj 

调节 项 为 XjZy (见方 程 (3)), 调节 效应 的 大 小 由 
归 系 数 y30 表示 ， 如 果 回 归 系 数 y 显著 不 为 0， 则 
表示 调节 效应 显著 ,简单 斜率 为 yioty302y。 需 要 
说 明 的 是 , 方程 2) 中 的 oops mays Hoy 和 uy 表示 层 
2 的 残 差 , 说 明 截 距 Bo, 和 斜率 Bij. Boy 和 pj 在 个 
体 间 可 以 变化 ( 称 为 随机 截 距 和 随机 斜率 )， 反 映 
了 个 体 间 的 差异 。 实 际 研究 中 ,为 了 简化 模型 ， 研 
究 者 可 将 Bij. Boy 和 py 设 为 固定 效应 ( 即 令 uij uzy 
和 Ha 为 0)。 

方程 (D)~(3) 所 示 的 纵向 数据 的 调节 效应 分 析 
不 仅 操作 方便 ， 而 且 容 易 解释 , 方便 应 用 。 但 这 种 
方法 有 不 足 之 处 ， 即 层 1 变量 将 随时 间 变 化 的 
个 体内 (intra-individual) 部 分 X, - X., 和 不 随时 间 
变化 的 个 体 间 (inter-individual) 部 分 了.; 混在 一 起 ， 
同 理 Zy 将 随时 间 变 化 的 个 体内 部 分 Z, -2.; 和 不 
随时 间 变 化 的 个 体 间 部 分 Z., 混在 一 起 (方术 等 ， 
2018; Preacher et al., 2016)。 
12 无 混淆 的 多 层 调节 效应 分 析 
121 显 变量 的 多 层 调节 效应 分 析 

为 了 解决 前 述 混 消 问 题 , Preacher 等 (2016) 建 


议 用 无 混淆 的 多 层 模型 (unconflated multilevel 

model) 进 行 调节 效应 分 析 ， 即 将 层 1 自 变量 Xy 
Zy 按 组 均值 E., WZ, 中 心 化 后 进行 建 模 ， 同 时 
层 -1: Yj = Poj + Bj (Xy-X.))+ BZ, -2.)) 
+f (X, -X,\(Z, —Z,)) + E; (4) 


层 -2: Poj = Yoo + yorX.; + Yo2Z.; + 
YosX.j2.; + Ho; (5a) 
Bij = V10 + 711Z. + A (5b) 
Paj =Y tho; 2) =o + 44; (Sc) 
将 方程 (5a)~(5c) 代 入 方程 (4) 得 : 
Y; = o + Nol Xy -X,,) + 7X; +Yo9(Z,y -Z,;) + 
Yor. + ¥59( Xj — XZ y —Z.j) + Yo3X.jZ.; + 
Ni (Xy -XZ + Ej + Moy +4; (Xy— Xj) + 
Li(Zy -2.)) + 1 (X; -X.;XZ; -Z.;) (6) 
其 中 ,系数 yo 只 表示 变量 Xy 随时 间 变 化 的 个 体 
内 部 分 x, -部 ， 的 效应 ,系数 yo 只 表示 变量 Xy 
不 随时 间 变 化 的 个 体 间 部 分 X., 的 效应 ， 从 而 区 
分 了 变量 Xy 的 个 体内 和 个 体 间 效应 (方程 (6) 与 方 
程 (3) 对 比 )。 同 理 ， 系 数 yo 只 代表 变量 Zy 随时 间 
变化 的 个 体内 部 分 Zy- Z. 的 效应 ， 系 数 yoo 只 表 
示 变 量 Zr 不 随时 间 变 化 的 个 体 间 部 分 Z,; 的 效应 ， 
从 而 区 分 了 变量 Zy 的 个 体内 和 个 体 间 效应 (方程 
(6) 与 方程 3) 对 比 )。 无 混 涌 ? 在 本 文中 的 含义 是 指 ， 
可 以 区 分 变量 的 个 体内 和 个 体 间 效应 。 
纵向 数据 的 调节 效应 也 分 为 三 个 部 分 (图 
1(b)) 。 第 一 部 分 是 随时 间 变 化 的 个 体内 部 分 
x, —¥,, W Z,- Z. 的 调节 效应 (以 下 称 纯 个 体内 
x 个 体内 (简写 为 WxW) 调 节 效 应 )， 调节 项 为 
X; -X.Z -2.))， 调 节 效 应 的 大 小 由 回归 系 
数 yo 表示 ， 如 果 回 归 系 数 y30 显著 不 为 0， 则 表 
示 纯 WxW 调节 效应 显著 ,简单 斜率 为 
yioty30(Z, -Z.,)。 第 二 部 分 是 不 随时 间 变 化 的 个 
EERS Æ., AZ, 的 调节 效应 (以 下 称 纯 个 体 间 
x 个 体 间 (简写 为 BxB) 调 节 效 应 )， 调节 项 为 
ZZ ,调节 效应 的 大 小 由 回归 系数 yn 表示 ,如 
果 回 归 系 数 yo3 显著 不 为 0， 则 表示 纯 BxB 调节 
效应 显著 ,简单 斜率 为 ortho 五 ,。 第 三 部 分 是 随 
时 间 变 化 的 个 体内 部 分 X, - Æ., 和 不 随时 间 变化 
的 个 体 间 部 分 Z.; 的 跨 层次 调节 效应 (以 下 称 个 体 
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间 x 个 体内 (简写 为 BxW) 调 节 效 应 )， 调 节 项 为 
(Xj —X.)Z.,; ， 调 节 效 应 的 大 小 由 回归 系数 yu 表 
示 ， 如 果 回 归 系 数 y 显著 不 为 0， 则 表示 跨 层 次 
调节 效应 显著 ,简单 斜率 为 yioty11 2.;° 
需要 说 明 的 是 , 方程 5(c) 的 pp 和 py 设 为 随机 
斜率 便于 研究 者 根据 实际 情况 添加 解释 变量 。 例 
wm, 将 Pa; =7%0+ bj 改 为 pi = Ya + YX. + Haj» 
则 增加 一 个 跨 层 次 调节 效应 ， 调 节 项 为 (Z, - 
2Z.;)X.;， 如 果 回 归 系 数 yo 显著 不 为 0， 则 表示 调 
节 变量 亏 ， 对 (Z, - 丈 ,) 一 关系 的 跨 层次 调节 效 


pee 另外 ， 由 于 无 混淆 的 多 层 模 型 将 分 为 
X,-X., MX, 两 个 变量 ,， 则 可 以 将 X.; 当成 调 


E ;对 (Xj 一 了 .,) 一 Yj 关系 的 调节 
效应 ， 此 时 并 没有 一 个 新 的 调节 变量 ZLIB, 
刘 红 云 , 2020)。 
1.2.2” 潜 变量 的 多 层 调节 效应 分 析 
基于 多 层 模 型 的 纵向 数据 的 调节 效应 检验 中 
(方程 (4)~(6)), 需要 使 用 不 随时 间 变 化 的 个 体 间 
部 分 X., 和 2Z.;。 多 层 模 型 设 定 变 量 为 显 变 量 ， 即 
默认 1 次 重复 测量 的 Xj 和 妈 2 的 均值 就 等 于 XX,; 和 


Level 2 


Levell 1 


图 2 基于 多 层 结构 方程 模型 的 纵向 数据 的 调节 模型 


进行 纯 WxW 的 层 1 调节 效应 分 析 ， 其 中 纯 WxW 
调节 项 为 RyRaj, 如 果 回 归 系 数 y 显著 ， 就 表示 
纯 WxW 调节 效应 显著 ( 见 网 络 版 附录 程序 )。 第 三 ， 
BxW 的 跨 层次 调节 效应 分 析 可 用 两 种 方式 进行 
分 析 。 一 种 是 随机 系数 预测 法 (Random Coefficient 
Prediction), 将 方程 (5b) 变 为 Bi; = 70 +U +4; > 
做 Bij 对 Uj 的 回归 ， 如果 回 归 系 数 yl 显著 ,就 表 
示 跨 层次 调节 效应 显著 ( 见 网 络 版 附录 程序 方 杰 
等 , 2018; Preacher et al., 2016)。 男 一 种 是 潜 调 节 
结构 方程 法 ， 利用 Rijs Rejs Uy 和 调节 项 RU, 
进行 BxW 的 跨 层 次 调节 效应 分 析 。 就 跨 层次 调节 

节 


T 


Z.;, BA N 次 重复 测量 中 抽取 部 分 重复 测量 结 
FR(tSN) KK X., AZ, 会 产生 抽样 误差 ( 方 杰 ， 
HARE, 2011; 方 杰 等 2014, 2018)。 比 较 好 的 解 
决 方法 是 将 并, 和 2.; 看 成 是 无 法 直接 观测 的 潜 
变量 ， 使 用 多 层 结构 方程 模型 进行 分 析 ， 能 较 好 
控制 X., MZ., ;的 抽样 误差 ( 方 杰 等 ，2018; Preacher 
et al., 216. 

结构 方程 模型 的 纵向 数据 的 调节 效应 分 
ris 2 2 所 示 。 第 一 , 将 每 个 层 1 测量 的 观测 变 
量 都 分 解 为 个 体内 部 分 和 个 体 间 部 分 两 个 潜 变 量 
CAS, EBA EE, 2011; FA, 2014, 2018; 
Preacher et al., 2016)。 例 如 ,观测 变量 Xy 分 解 为 
个 体内 部 分 潜 变 量 Rey 和 个 体 间 部 分 潜 变 量 Uj。 
同 理 , Zj 分 解 为 两 个 潜 变 量 Roy 和 Uy, 7 分 解 为 两 
个 潜 变 量 Ry 和 U,。 第 二 ,使 用 潜 调 节 结 构 方程 
LMS, 


AS 


法 (Latent Moderate Structural Equations, 
Klein & Moosbrugger, 2000), 根据 

Uy = Yo t+ VU y + YoU y + YoU yy + Ho; (7) 
进行 纯 BxB 的 层 2 调节 效应 分 析 ， 其 中 纯 BxB 调 
节 项 为 UyU;y， 如果 回归 系数 yo 显著 ， 就 表示 纯 
BxB 调节 效应 显著 。 同 理 , 根据 

Ry = Boj + Nowy + yr0Rzy +y30RwyRzy + Ey (8) 


效应 分 析 而 言 ,随机 系数 预测 法 的 计算 量 比 潜 调 
结构 方程 法 要 少 ( 方 杰 等 ，2018; Preacher et al., 
2016)， 所 以 建议 选用 随机 系数 预测 法 ; 进行 路 层次 
调节 效应 分 析 。 


2 基于 交叉 滞后 模型 的 纵向 数据 的 调 
节 效 应 分 析 


2.1 ”传统 的 交叉 滞后 的 调节 效应 分 析 
传统 的 交叉 滞后 模型 的 调节 效应 分 析 ( 图 3) 
可 表示 为 
Xi = ByX i(t-1) +7 xX (t- 1) + Cx Xie- 1) Yiu 1) + Ex (9) 
= Bua) + WX iq 1) + Cy Xili) + Evu (10) 


其 中 , 下 标 i 表示 个 体 ， 系数 px 和 py 表示 自 回归 
效应 (autoregressive effect， 即 每 个 变量 前 后 时 间 
点 的 效应 )， 系数 yw 和 条 表 示 灌 后 效应 (lag effect, 
即 每 个 变量 对 前 因 变 量 的 回归 )。 调 节 项 是 
Xa- WRR Cr MER) Cy 显著 不 为 0， 则 
表示 调节 效应 显著 。 值 得 注意 的 是 ， 自 回归 系数 、 
滞后 系数 和 调节 系数 的 下 标 都 没有 时 间 t, 表示 
月 回归 效应 、 滞 后 效应 和 调节 效应 在 不 同时 间 点 
都 具有 稳定 性 (目的 是 简化 模型 )。exr 和 sz 表示 误 
差 项 ， 服 从 正 态 分 布 。 
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图 3 传统 的 交叉 滞后 调节 模型 


交叉 滞后 模型 的 调节 效应 至 少 有 两 个 意义 
(以 cx 为 例 进行 说 明 )。 第 一 ,所 调节 了 滞后 效 
应 Yia-) Xit o 第 二 ， Yio al T 自 回 归 效 应 
Xia- 7X, (Ozkok et al., in press), HAZ, 26M HA 
模型 可 以 考察 具有 时 间 先 后 关系 的 变量 间 的 调节 
效应 ， 这 是 交叉 滞后 模型 的 一 大 特色 。 值 得 说 明 
的 是 , 在 截面 数据 的 调节 研究 中 ,调节 变量 Z 是 
一 个 不 同 于 自 变 量 X 和 因 变 量 了 的 变量 。 但 是 , 在 
纵向 数据 的 调节 研究 中 ， 就 不 一 定 是 这 样 了 。 例 
如 , 在 交叉 滞后 模型 中 ， 自 变量 筷 和 因 变 量 了 变 成 
了 四 个 变量 CC Xen. Ya 站 D)， 就 可 以 不 需 
要 一 个 新 的 调节 变量 Z 了 (公式 (9) 和 (10))， 这 是 纵 
向 数据 的 调节 研究 有 别 于 截面 数据 的 一 个 特点 。 
2.2 无 混淆 的 交叉 滞后 的 调节 效应 分 析 
传统 的 交叉 滞后 模型 的 纵向 数据 的 调节 分 析 
方法 (见方 程 (9) 和 (10)) 同 样 存 在 不 足 ， 即 变量 (如 
Xa M Z) 也 将 随时 间 变 化 的 个 体内 部 分 和 不 随时 
间 变 化 的 个 体 间 部 分 混在 一 起 。 为 了 解决 这 一 不 
JE, Hamaker 等 (2015) 建 议 将 中 心 化 后 的 观测 分 数 
Xi DRA RE Xp; 和 Xm (图 和 ， 即 
Xu = Xa + Kyi (11) 
其 中 , Xg: 表示 不 随时 间 变 化 (下 标 没 有 时 间 EY 
本 间 差异 ， 即 随机 截 距 ; Yim 表示 随时 间 变 化 的 个 
林内 差异 ， 从 而 将 变量 总, 的 个 体内 和 个 体 间 部 分 
区 分 开 。 同 理 ,将 中 心 化 后 的 观测 分 数 歹 分 解 为 
两 个 潜 变 量 Yo; 和 Ym (Al 4)， 即 
Yi =YB; + Yr (12) 


从 而 将 变量 二 ,的 个 体内 和 个 体 间 部 分 区 分 开 。 此 
时 , 方程 (9) 和 (10) 的 调节 项 Xi Yee : 
Xia pho n= (Xa + Xap) * 

(Ya: + Yria-)) = 

X gia + X gima- + 


a3) 


YgiX mia- + Xi viel) 


可 以 看 到 X41)Yio-y 实 际 上 包含 了 3 种 调节 效应 
(Ozkok et al., in press)。 第 一 种 是 纯 WxW 调节 效 
应 ， 即 随时 间 变 化 的 个 体内 调节 效应 
Xme-Ymen; 第 二 种 是 纯 BxB 调节 效应 ， 即 不 随 
时 间 变 化 的 个 体 间 调节 效应 X57s， 第 三 种 是 
BxW 跨 层 次 调节 效应 Ya Xie Al XaiYma-1)。 接 下 
来 将 分 别 介绍 这 三 种 调节 效应 分 析 , 然后 再 介绍 
三 种 调节 模型 的 整合 分 析 。 
2.2.1 ” 纯 个 体内 的 潜 调 节 效 应 分 析 

纯 个 体内 的 潜 调 节 效 应 分 析 ( 图 4) 可 表示 为 
方程 (14) 和 (15) 

Xwi = BxXma-y + YxYma- + 
CyX ma) Ymca + Exi 
Ya = PrY wig ay + YrXwia-) + 
Cy X ma Yi + Evin 

其 中 , HEIA Xmen Ymen WRR Cr PER) 
Cy 显著 ,表示 调节 效应 显著 (Ozkok et al., in 
press). Ozkok 等 (in press) 用 模拟 研究 证 明 ， 纯 个 
体内 的 潜 调 节 效 应 分 析 ( 方 程 (14) 和 (15)) 优 于 传统 
的 交叉 滞后 的 调节 效应 分 析 ( 方 程 (9) 和 (10))， 表 
现 为 更 小 的 偏差 估计 和 更 大 的 收敛 率 。Ozkok 等 
还 用 模拟 研究 发 现 ， 纯 个 体内 的 潜 调 节 效 应 分 析 
优 于 纯 个 体内 的 显 变量 调节 效应 分 析 ( 即 方程 (11) 
PÉ Xman = Xop tB X, Xai = Xu; PIM, 


方程 (12) 中 的 Yma- = Yio Ya, Ypi = Fo Jo 


a4 


a5) 


图 4 ai Wx W 的 交叉 滞后 调节 模型 


(改编 自 Ozkok et al., in press) 


2.22 ” 跨 层 次 的 潜 调 节 效 应 分 析 
跨 层 次 的 淤 调节 效应 分 析 ( 图 5) 可 表示 为 方 
程 (16) 和 (17) 
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Xyu = BxXwia-y + YxYwia- + CixX piYwicn) + 

Coy XBpX win + CaxYBpX wia-) + 

Cax YBpiY mia- + Exi (16) 

Yi = PY wey + YrXmia- + ClyX aX wie—y + 

CyyX eX iat) + Coy¥ aX wig ay + 

Cay Ypi Yma- + Evin (17) 
其 中 ,调节 项 为 XpiYmia 1)、XBAXmio 1)、YBiX ma1) 
和 YoY piety, 如果 Cix、 Cox. Cax、 Cax、 Crys Coys 
Gy, Car 中 的 某 个 系数 显著 ,就 表示 与 该 系数 相 
应 的 调节 效应 显著 。 与 方程 (14) 和 (15) 相 比 ， 跨 层 
次 的 潜 调 节 效 应 模型 (方程 (16) 和 (17)) 仅 仅 是 调节 
项 变 多 了 。Ozkok 等 (in press) 用 模拟 研究 证 明 ， 跨 
层次 的 潜 调 节 效 应 分 析 ( 方 程 (17) 仅 保留 调节 项 
YsiXwmia1) 也 优 于 跨 层次 的 显 变量 调节 效应 分 析 
( 即 Ygi = Ya ` Xma- = Xey X), 表现 为 Coy HY 
偏差 估计 更 小 和 更 大 的 收敛 率 。 


< 
AES 
D 


图 5 BxW 的 交叉 滞后 调节 模型 
(改编 自 Ozkok et al., in press) 
2.2.3 ” 纯 个 体 间 的 潜 调节 效应 分 析 
纯 个 体 间 的 潜 调 节 效 应 分 析 可 表示 为 方程 
(18) 和 (19) 
AP = DyZ pi +YxYpi + CyZeiY ai + Ex (18) 
Yp; = DyZp; + Vy X pi + CyZ pi X gi tEn (19) 


其 中 , Zp; 是 一 个 新 的 个 体 间 变 量 (如 性 别 )， 调 节 
项 是 ZgiYpi 和 ZpiXp;， 如 果 Cy M(BOCy MH, 表示 
调节 效应 显著 (Ozkok et al., in press)。 
2.2.4 ”三 种 潜 调 节 效 应 的 整合 分 析 

如 果 一 个 研究 只 有 两 个 变量 ,和 Yi, 可 以 将 
两 种 潜 调 节 效 应 整合 在 一 起 进行 分 析 。 具 体 地 ， 
由 方程 (14) 和 (15) 进 行 个 体内 的 潜 调 节 效 应 分 析 ， 
方程 (16) 和 (17) 进 行 跨 层 次 的 潜 调 节 效 应 分 析 即 
可 , 无 法 进行 个 体 间 的 潜 调 节 效 应 分 析 。 


如 果 自 变量 Xy 和 因 变 量 ,随时 间 变 化 (纵向 
数据 )， 调 节 变 量 Z5; 不 随时 间 变 化 (截面 数据 )， 就 
可 以 将 三 种 潜 调 节 效 应 整合 在 一 起 进行 分 析 。 具 
体 地 ， 由 方程 (14) 和 (15) 进 行 个 体内 的 潜 调 节 效 应 分 
Br, 方程 (18)~(20) 进 行 个 体 间 的 潜 调 节 效 应 分 析 ( 如 
果 调 节 变 量 是 二 分 类 别 变量 , 则 去 掉 方程 (20))， 

Zp = DY pi + YzX pi + CXBpiYBi + Ez (20) 
跨 层次 的 潜 调 节 效 应 分 析 则 需要 在 方程 (16) 和 (17) 
中 都 各 自 加 上 两 个 调节 项 ZpiYwo_1) 和 Ze X wie El 
可 (Ozkok et al., in press)。Ozkok 等 人 用 一 个 示例 
演示 如 何 将 潜 调 节 结 构 方 程 法 和 Bayesian 法 相 结 
合 ， 对 不 随时 间 变 化 的 调节 变量 (性 别 ) 做 三 种 潜 
调节 效应 的 整合 分 析 。 

如 果 所 有 变量 都 随时 间 变 化 (如 Xi. Vir FA Zi), 
同样 也 可 以 将 三 种 潜 调 节 效 应 整合 在 一 起 进行 分 
析 。 第 一 , 方程 (14) 和 (15) 都 各 自 增加 Zmes 
Xmo_DZmD 和 Yo_wZmen) 三 项 ， 再 加 方程 (21)， 
即 可 进行 个 体内 的 潜 调 节 效 应 分 析 。 

Lwin = PzZLwia ny) + YzxX mi + Yzy ľwia- + 

CyzX wia-)Zmia-) + Cry X moa pyYwia-n + 

Cyz Yi Z wie + Ezi (21) 
第 二 ， 由 方程 (18)~(20) 进 行 个 体 间 的 潜 调 节 效 应 
分 析 。 第 三 ， 方程 (16) 和 (17) 中 都 各 自 增加 Zwict-1) 
ZBiYma-1) > ZeXma-1) > ZeiZwi-vy > YeiZwic-t) > 
XBL wi(t-1)s 再 加 上 方程 (22)， 即 可 进行 跨 层次 的 潜 
调节 效应 分 析 。 

Zwi = B22 wie) + Vax wia- + Yzy wia- + 
C2 Va X wia- + CozYBpiYwic nt) + 


C32YgiZwia-) + CazX rimia- + 

CszX pi X mia- + CozX g:Zwia-) + 

C57 Z ayia) + CezZ wi X wia- + 

CozZgiZwia-1) + Ezi (22) 
需要 说 明 的 是 ， 三 种 潜 调 节 效 应 的 整合 分 析 会 有 
很 多 个 调节 项 ， 这 无 疑 增加 了 数据 不 收敛 的 可 能 ， 
也 大 大 增加 了 运行 时 间 。 在 实践 中 ,建议 研究 者 
先 从 理论 上 排除 一 些 没有 充分 理论 依据 的 调节 项 ， 
从 而 减少 计算 时 间 。 


3 基于 潜 变 量 增 长 模型 的 纵向 数据 的 
调节 效应 分 析 
如 果 自 变量 、 因 变量 和 调节 变量 都 随时 间 变 
化 (纵向 数据 )， 则 潜 变 量 增长 模型 的 纵向 数据 的 
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图 6 基于 潜 变 量 增长 模型 的 纵向 数据 的 调节 模型 (改编 自 Wen et al., 2014) 


调节 效应 分 析 包括 四 个 潜 变 量 增长 模型 ， 分 别 是 
自 变 量 的 潜 变 量 增长 模型 、 调 节 变 量 的 潜 变 量 增 
长 模型 、 因 变量 的 潜 变 量 增长 模型 、 和 调节 项 的 
潜 变 量 增长 模型 (图 6), 分 析 的 关键 是 如 何 构 建 调 
节 项 的 潜 变 量 增长 模型 , Wen 等 (2014) 采 用 乘积 指 
标 构 建 了 调节 项 的 初始 水 平 ,用 差 积 指标 构建 了 
调节 项 的 变化 率 。 具 体 地 ， 自 变量 的 潜 变 量 线 性 


Hee RID SS 


增长 模型 可 以 表示 为 
Xi=ar+O (23a) 
X,=axy +Sy +ô, (23b) 
X3=ax +25% +06; (23c) 
X=ay +3Sy +ô, (23d) 


其 中 , ax 和 Sx 都 是 潜 变 量 , 4 AREE X WY 
水 平 ( 截 距 ) 和 变化 率 。 同 理 , 调节 变 量 的 潜 变 量 
线性 增长 模型 可 以 表示 为 


Zi=az +ô; (24a) 
Zy =a, + Sy +ô (24b) 
Z,=az +2S7 +ô (24c) 
Z4=az +357 +ô (24d) 
将 方程 (23a) 乘 以 (24a)， 整 理 得 : 
XZ, =aydz +70, +axos + 6,05 (25) 


至 此 ,初始 水 平 的 调节 效应 模型 (图 6)， 即 az 对 ax 
和 ay 关系 的 调节 效应 可 表示 为 

ay=b+ day +eaz + fa yaz + E,y (26) 
其 中 , 调节 项 为 潜 变 量 axaz X\Z, 是 ayaz 的 指标 ， 
由 于 是 单 指标 潜 变 量 ， 因 此 固定 其 负 蓓 为 1。 如 果 
系数 了 显著, 表示 初始 水 平 的 调节 效应 显著 ( 温 忠 


BE, XLA, 2020; Wen et al., 2014)。 
将 方程 (23b) 减 去 (23a) 、 方 程 (24b) 减 去 (24a)， 
得 到 : 


(X, — X)) =Sy + (6, — 6) (27a) 
(Z, —Z,) = Sz + (6% — 5s) (27b) 
将 方程 (27a) 乘 以 (27b)， 整 理 得 ; 
(X2 -X22 五 )=9x9z+Sxr(05 一 05) 二 
Sz(0 — 61) + (8> -Â s —6;) (28) 
同 理 ， 由 方程 (23d)、(23c) 和 方程 (24d)、(24c), 得 到 : 
(Xa—X3)(Z4 —Z3) =SySz +Sy(ð; - 0) + 
S7(, — 63) + (8, -6% — 67) 


至 此 ,变化 率 的 调节 效应 模型 (图 6), 即 Sz 对 Sy 
和 Sy 关系 的 调节 效应 可 表示 为 
Sy=g +hSy +kSz + wSySz + Esy (30) 
其 中 , 调节 项 为 潜 变 量 SySz,(X, —X{)(Z, -Z,) Fil 
(X,-X3)(Z,-Z,) 是 Sx5z 的 指标 (图 6)， 如 果 系 
数 w 显著 ,表示 变化 率 的 调节 效应 显著 ( 温 忠 甩 ， 
刘 红 云 , 2020; Wen et al., 2014)。 
除了 上 述 需 要 构建 乘积 项 的 乘积 指标 法 ,还 
可 以 使 用 无 需 构 建 乘积 项 的 潜 调 节 结 构 方 程 法 进 
行 基 于 潜 变 量 增长 模型 的 纵向 数据 的 调节 效应 分 
析 ( 详 见 第 6 部 分 的 示例 )。 潜 调节 结构 方程 法 的 一 
个 优势 是 , 除了 初始 水 平 (调节 项 axaz) 和 变化 率 
(调节 项 Sx5z) 的 调节 效应 外 ， 潜 调节 结构 方程 法 
还 可 以 便利 地 产生 其 它 调节 项 (如 axSx、axSz 和 
azS)， 灵 活 地 进行 其 它 调 节 效 应 分 析 。 甚 至 当 调 
节 变 量 Z 是 不 随时 间 变 化 时 , 潜 调 节 结 构 方程 法 
也 可 以 便利 地 产生 调节 项 (如 axZ、SxZ)， 进行 调节 
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效应 分 析 。 男 外 ， 由 于 潜 变 量 增长 模型 将 自 变量 外 
分 为 ax 和 Sx 两 个 潜 变 量 ， 从 理论 上 说 ， 则 可 以 将 
Sy 当成 调节 变量 ,利用 潜 调 节 结 构 方 程 法 分 析 Sx 
对 az 一 Sr 关系 的 调节 效应 ， 并 不 一 定 需要 一 个 新 
的 调节 变量 Z。 


4 纵向 数据 的 调节 模型 的 综合 比较 


多 层 模 型 (MLM)、 多 层 结构 方程 模型 (MSEMD)、 
交叉 滞后 模型 (CLM) 和 潜 变 量 增长 模型 LGMD) 的 
调节 模型 比较 如 表 1 所 示 。 纵 向 数据 的 调节 效应 
分 析 的 一 大 特点 是 调节 效应 的 分 解 。 具 体 地 ， 多 
层 模 型 、 多 层 结构 方程 模型 和 交叉 滞后 模型 都 将 
变量 区 分 为 随时 间 变 化 的 个 体内 部 分 (W) 和 不 随 
时 间 变 化 的 个 体 间 部 分 (B)， 从 而 将 纵向 数据 的 调 
节 效 应 分 解 为 三 部 分 ， 即 纯 WxW 调节 、 纯 BxB 


方便 实现 (Muthén & Muthén, 1998-2017). (44, 
Mplus 默认 在 潜 调 节 结 构 方 程 法 的 调节 效应 估计 
H, 采用 极 大 似 然 估计 (Muthén & Muthén, 1998- 
2017)。 当 自 变量 和 调节 变量 都 是 潜 变 量 时 , 极 大 
似 然 估计 必须 用 到 数值 积分 mumerical integration), 
这 将 导致 计算 量 很 大 (Asparouhov & Muthén, 
2021)。 有 研究 者 建议 ， 当 自 变 量 和 调节 变量 都 是 
潜 变 量 时 , 将 潜 调 节 结 构 方 程 法 和 Bayesian 法 相 
结合 (Mplus8.3 及 以 后 版 本 方 可 实现 )， 得 到 调节 
效应 的 Bayesian 佑 计 (Asparouhov & Muthén, 2021; 
Ozkok et al., in press). Asparouhov 和 Muthén(2021) 
的 模拟 研究 表明 ,在 渤 调 节 结 构 方程 法 的 调节 效 
应 分 析 中 ， 相 比 Preacher 等 (2016) 使 用 的 极 大 似 
TMA, Bayesian 估计 会 更 快 、 更 准确 (表现 
为 更 小 的 绝对 偏差 (absolute bias) AH HE UF Ay KX lal BE 


调节 和 BxW 路 层次 调节 。 潜 变量 增长 模型 将 纵向 
数据 的 调节 效应 分 解 为 初始 水 平 的 调节 效应 和 变 
化 率 的 调节 效应 。 

纵向 数据 的 调节 分 析 的 另 一 个 特点 是 进行 潜 
变量 的 调节 效应 分 析 , 而 且 大 都 推荐 使 用 潜 调 节 
结构 方程 (LMS) 法 进行 潜 变 量 的 调节 效应 分 析 
(Asparouhov & Muthén, 2021; Ozkok et al., in press; 
Preacher et al., 2016; Wen et al., 2014)。 具 体 地 ， 
Preacher 等 (2016) 明 确 指出 ， 只 有 潜 调 节 结 构 方 程 
法 能 适用 于 1 x (1 一 1) 调 节 模 型 的 WxW 调节 效应 
分 析 ( 即 方程 (8) 的 RRzy)o Wen 等 (2014) 比较 了 乘 
积 指 标 法 和 潜 调 节 结 构 方 程 法 在 潜 变 量 增长 模型 
的 变化 率 的 调节 效应 检验 中 的 表现 。 结 果 发 现 ， 
在 正 态 条 件 下 , 尤其 是 在 小 样本 时 ， 湾 调节 结构 
方程 法 比 乘积 指标 法 在 参数 估计 时 产生 更 小 的 偏 
差 。 但 是 在 数据 非 正 态 条 件 下 ,乘积 指标 法 比 潜 
调节 结构 方程 法 在 参数 估计 时 产生 更 小 的 偏差 。 
因此 ， 当 数据 正 态 分 布 条 件 下 ,推荐 使 用 潜 调 节 
结构 方程 法 ; 数据 非 正 态 条 件 下 ,推荐 使 用 乘积 
指标 法 。 

目前 ， 湾 调节 结构 方程 法 可 通过 Mplus 软件 


盖 率 (interval coverage rate) )、 且 收敛 率 更 高 。 

第 三 ,收集 纵向 数据 进行 分 析 的 一 个 重要 月 
的 就 是 为 了 进行 历时 性 因果 推断 ( 温 忠 腾 ，2017)， 
如 果 要 进行 历时 性 因果 推断 ， 在 本 文 介绍 的 几 类 
纵向 模型 中 ， 只 有 交叉 滞后 模型 及 其 变 体 可 以 实 
现 ,， 因 为 只 有 交叉 滞后 模型 体现 了 历时 性 因果 推 
测 所 必须 的 测量 时 间 点 先后 关系 ,同时 考虑 了 自 
回归 效应 ( 方 杰 等 , 2021)。 其 他 三 种 模型 在 建 模 
过 程 中 , 没有 考虑 变量 之 间 影 响 的 先后 顺序 ， 其 
调节 效应 分 析 本 质 上 与 截面 数据 的 调节 分 析 无 异 ， 
要 得 出 因果 推论 还 要 有 统计 以 外 的 理据 ( 温 忠 麟 ， 
2017; 方 杰 等 , 2021)。 

另外 ， 如 果 是 显 变量 的 调节 效应 检验 或 想 耗 
时 少 ， 建 议 选择 多 层 模型 ; 如 果 是 测量 次 数 只 有 
两 次 ， 建 议 选择 交叉 滞后 模型 ; 如 果 是 考察 初始 
水 平和 (或 ) 变 化 率 的 调节 效应 检验 ,建议 选择 洪 


变量 增长 模型 。 
5 纵向 数据 的 调节 效应 分 析 流程 


看 对 一 个 纵向 数据 的 调节 效应 分 析 任 务 ， 研 
究 者 应 当 如 何 进行 呢 ? 根据 前 面 的 讨论 ， 我 们 总 


表 1 纵向 数据 的 调节 模型 比较 


用 LMS 法 耗 时 


模型 潜 变 量 调节 调节 效应 的 分 解 测量 次 数 ”做 历时 性 因果 推断 

MLM x WxW、BxB 和 BxW 调节 三 3 x x 短 
MSEM v WxW、BxB 和 BxW 调节 =3 x V 长 
CLM vV WxW、BxB 和 BxW 调节 三 2 Vv vV K 
LGM v 初始 水 平 调节 、 变 化 率 调节 等 =3 x V k 


注 :“W” 表 示 可 以 考察 ,“x” 表 示 无 法 考察 。 


"202303.09505v1 


chinaXiv 


ChinaxX iv ERAT 


2468 心理 科学 进展 第 30 卷 


结 出 一 套 纵 向 数据 的 调节 效应 分 析 流 程 ( 见 图 7): 

(1) 需要 进行 历时 性 因果 推论 吗 ? 如 果 是 ， 
则 用 交叉 滞后 模型 进行 调节 效应 分 析 。 否 则 ， 进 
人 步骤 2。 

(2) 将 纵向 数据 当成 多 层 数 据 吗 ?如果 是 ， 则 
用 多 层 结构 方程 模型 进行 调节 效应 分 析 。 多 层 结 
构 方 程 模型 和 多 层 模型 更 适合 于 描述 个 体 间 差 
异 。 否 则 , 使 用 潜 变 量 增长 模型 进行 调节 效应 分 


(2013) 的 研究。 采用 Mplus 8.3 软件 进行 纵向 数据 
的 调节 效应 分 析 (Mplus 程序 见 网 络 版 附录 )。 
6.1 基于 交叉 滞后 的 纵向 数据 的 调节 效应 分 析 
如 果 需 要 进行 历时 性 因果 推论 ， 则 使 用 交叉 
滞后 模型 进行 纵向 数据 的 调节 效应 分 析 。 为 了 简 
化 模型 ,本 例 使 用 家 庭 冲突 (XY) 和 青少年 被 欺负 (7) 
两 个 变量 , 将 Bayesian 法 和 潜 调 节 结 构 方程 法 相 结 
合 进行 纯 WxW 的 无 混淆 的 交叉 滞后 模型 的 潜 调 节 


析 。 只 有 潜 变 量 增长 模型 可 以 同时 研究 某 些 变量 
对 变化 的 影响 (将 变化 作为 因 变 量 ,， 如 Sy) 和 变化 
如 何 影 响 其 他 变量 ( 即 把 变化 作为 自 变量 , 如 Sy) 
(iL BE, 刘 红 云 , 2020)。 

(3) 多 层 结构 方程 模型 收敛 吗 ? 如果 是 , 报告 
多 层 结 构 方程 模型 的 调节 效应 分 析 结 果 。 否 则 ， 
使 用 多 层 模型 进行 调节 效应 分 析 。 多 层 结构 方程 模 
型 和 多 层 模 型 的 主要 差别 在 于 计算 组 均值 时 是 否 
考虑 抽样 误差 ( 见 1.2.2 部 分 ), 但 多 层 结 构 模 型 的 
收敛 难度 更 大 ， 因 此 优先 使 用 考虑 了 抽样 误差 的 
多 层 结 构 方 程 模型 ， 收 敛 不 了 再 考虑 多 层 模 型 。 


需要 进行 历时 性 因果 推论 吗 ? 


图 7 纵向 数据 的 调节 效应 分 析 流 程 图 


6 纵向 数据 的 调节 效应 分 析 示 例 


下 面 用 一 个 实际 例子 演示 如 何 用 图 7 的 流程 
进行 纵向 数据 的 调节 分 析 。 本 例 的 变量 有 家 庭 冲 
突 (V)、 青 少年 被 欺负 (7) 和 烟 酒 使 用 ， 变量 及 其 数 
据 ( 每 间隔 半年 重复 测量 1 次 , 共 重 复 测量 4 次 ， 
FE 1132 名 11~15 岁 的 美国 青少年 ) 均 来 自 Little 


效应 分 析 ( 方 程 (14) 和 (15)， 有 3 个 调节 项 XlYmi、 
XwaYm 和 XmaYm3)。 分 析 前 ， 自 变量 和 因 变 量 按 
总 均值 进行 了 中 心 化 。 

第 一 ，Bayesian 法 是 以 马尔 科 夫 链 (Markov 
Chain) 收 敛 为 前 提 的 ,常用 轨迹 图 (trace plot， 如 
果 多 条 马尔 科 夫 链 的 轨迹 高 度 重 亚 表 明 马 尔 科 夫 
链 收 敛 ) 和 PSR(Potential Scale Reduction， 如 果 
PSR<1.1 (或 更 严格 的 标准 1.05) 表 明 马 尔 科 夫 链 
收敛 ) 来 检验 马尔 科 夫 链 是 否 收敛 , 更 详细 的 介绍 
见 综述 方 杰 和 张 敏 强 等 (2011)。 轨 迹 图 和 PSR 结 
果 (TECHS8 报告 第 7700 次 迭代 时 PSR = 1.049) 表 
明 马 尔 科 夫 链 收 敛 。 

第 二 , 纯 WxW 的 潜 调 节 效 应 分 析 结 果 显 示 ， 
Xma- Xm W A EIIZ By = 0.482, SE = 0.025, 
95% 可 靠 区 间 [0.431, 0.532] 不 包含 0， 自 回归 效应 
By 显著 , 表明 随 着 时 间 的 推移 ， 家 庭 冲 突 逐 渐 加 
Eo Viney Xie Wi OD yy = 0.003, SE = 0.02, 
95% HY fk XE] [-0.037, 0.044] 0, 滞后 效应 yy 
不 显著 , 表明 1 时 刻 的 青少年 被 欺负 不 能 显著 
预测 后 一 时 刻 t 的 家 庭 冲突 。Xwa_DYma_D 王 Xi 
的 调节 效应 Cy = —0.053, SE = 0.026, 95% 可 靠 区 
间 [-0.105，-0.004] 不 包含 0, 调节 效应 Cx 显著 ， 
表明 Xp 调节 了 滞后 效应 Yme Xm 或 者 
Yma- T 自 回归 效应 Xm IX wits 

B=, Yme- Ym H A PIRU py = 0.286, 
SE = 0.028, 95% 可 靠 区 间 [0.228, 0.341] 不 包含 0， 
自 回 归 效 应 py 显著 , 表明 随 着 时 间 的 推移 ， 青 少 
年 被 欺负 逐渐 加 重 。 这 也 符合 实际 ,被 欺负 的 青 
少年 更 容易 继续 被 欺负 。 Xmio_w 一 Ym 的 滞后 效应 
yy = 0.07, SE = 0.028, 95% 可 靠 区 间 [0.016, 0.126] 
不 包含 0, 滞后 效应 yy 显著 , 表明 1 时 刻 的 家 庭 
冲突 显著 增加 了 后 一 时 刻 t 的 青少年 被 欺负 。 这 
可 解释 为 ,青少年 和 父母 冲突 越 多 , 会 让 其 他 人 
认为 ,这 个 青少年 缺少 安全 保护 , 会 更 容易 将 之 
作为 欺负 的 对 象 。Xywo-pYmio-y 一 Ym 的 调节 效应 
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Cy = -0.039，SE = 0.033，95% 可 靠 区 间 [-0.105， 
0.024] 包 含 0， 调 节 效 应 Cy 不 显著 。 
第 四 ,不 随时 间 变 化 的 个 体 间 关系 Xe Yri 
的 效应 为 0.564, SE = 0.058, 95% 可 靠 区 间 [0.459， 
0.681] 不 包含 0, Xp 一 Yi 的 效应 显著 ， 表明 家 庭 } 
突显 著 增加 了 青少年 被 欺负 的 风险 。 需 要 说 明和 
Fe, 本 例 的 纯 W x W 的 潜 调 节 效 应 检验 (有 3 个 】 
节 项 ) 耗 时 3 个 半 小 时 ， 因此， 跨 层次 的 潜 调 节 效 
应 检验 会 耗 时 更 长 。 
62 ”基于 多 层 结 构 方程 模型 的 纵向 数据 的 调节 
效应 分 析 


Kp 


= 


调节 结构 方程 法 相 结合 进行 基于 潜 变 量 增长 模型 
的 潜 调 节 效 应 分 析 ( 有 二 个 调节 项 axaz 和 SySz)o 
第 一 ,轨迹 图 和 PSR 结果 (TECH8 报告 第 
1900 次 迭代 时 PSR = 1.011) 表 明 马 尔 科 夫 链 收敛 。 
第 二 ,变化 率 的 调节 效应 检验 结果 发 现 ( 以 下 呈现 
的 是 标准 化 解 ), Sx 对 Sy 的 效应 显著 (h = 0.363, SE 
= 0.084, 95% 可 靠 区 间 [0.2, 0.528] 不 包含 0); Sz 对 
5y 的 效应 显著 (k = -0.21, SE = 0.105, 95% 可 靠 区 间 
[-0.418, 0.007] KH 0); 变化 率 的 调节 效应 显 
著 (w = 0.052, SE = 0.024, 95% 可 靠 区 间 [0.006, 0.098] 
不 包含 0)， 即 烟 酒 使 用 的 变化 率 正 向 调节 了 家 庭 


如 果 将 纵向 数据 当成 两 层 数 据 ， 则 使 用 多 层 
结构 方程 模型 进行 调节 效应 分 析 。 要 研究 的 是 


冲突 的 变化 率 与 青少年 被 欺负 的 变化 率 的 关系 。 
第 三 , 初始 水 平 的 调节 效应 检验 结果 发 现 ( 以 


内 
酒 使 用 ( 妃 对 家 庭 冲突 (XY) 和 青少年 被 欺负 (六) 的 调 
节 效 应 。 采 用 Bayesian 法 和 潜 调 节 结 构 方 程 法 、 
随机 系数 预测 法 相 结 合 进行 潜 调 节 效 应 分 析 。 其 
中 ,随机 系数 预测 法 仅 用 来 考察 跨 层次 的 潜 调 节 
效应 ， 其 余 潜 调节 效应 用 潜 调 节 结 构 方 程 法 进行 
考察 。 为 了 简化 模型 ,将 uy 和 jw3j 设 为 0。 
第 一 ,轨迹 图 和 PSR 结果 (TECH8 报告 第 900 
次 迭代 时 PSR = 1.044) 表 明 马 尔 科 夫 链 收敛 。 第 
二 ， 基 于 多 层 结构 方程 模型 的 纵向 数据 的 调节 效 
应 分 析 显 示 ， 纯 个 体内 ( 纯 WxW) 的 潜 调 节 效 应 
y30 = —0.191, SE = 0.077，95% 可 靠 区 间 [-0.35， 
一 0.047] 不 包含 0, y30 显著 不 为 0， 表明 纯 个 体内 的 
潜 调 节 效 应 显著 。 第 三 , 纯 个 体 间 ( 纯 B x B) 的 潜 
调节 效应 yos = 0.232, SE = 0.086, 95% 可 靠 区 间 
[0.063，0.394] 不 包含 0, y 显著 不 为 0， 表明 纯 个 


呈现 的 是 标准 化 解 )， ay 对 ay 的 效应 显著 (d = 
0.417, SE = 0.082, 95% 可 靠 区 间 [0.264, 0.581] 不 
包含 0); az 对 ay 的 效应 不 显著 (e = -0.117, SE = 
0.064, 95% 可 靠 区 间 [-0.251，0.001] 包 含 0); 初始 
水 平 的 调节 效应 显著 (f= -0.173, SE = 0.053, 95% 
可 靠 区 间 [-0.282，-0.074] 不 包含 0)， 即 烟 酒 使 用 
的 初始 水 平 调节 了 家 庭 冲突 的 初始 水 平 与 青少年 
被 欺负 的 初始 水 平 的 关系 。 


7 讨论 和 拓展 


相 比 截面 数据 ， 纵 向 数据 的 调节 效应 分 析 要 
复杂 的 多 。 对 于 自 变量 和 因 变 量 有 纵向 数据 的 调 
节 效 应 分 析 ， 可 根据 调节 变量 的 不 同类 型 分 为 三 
类 。 第 一 类 是 调节 变量 从 自 变量 或 因 变量 中 产生 。 
MEERA PIX, 对 (xX, -X.,) oY, 关系 的 调 


体 间 的 潜 调 节 效 应 显著 。 第 四 ， 跨 层次 的 调节 效 
DM yy = -0.113, SE = 0.156, 95% 可 靠 区 间 [-0.403， 
0.209) 0, y RAD, 表明 跨 层次 的 潜 调 节 效 
应 不 显著 。 其 它 参 数 估计 见 表 2。 
6.3 ”基于 潜 变 量 增长 模型 的 纵向 数据 的 调节 效 
应 分 析 
要 研究 的 是 烟 酒 使 用 (和 对 家 庭 冲 突 ( 妨 和 青 
少年 被 欺负 ( 尺 的 调节 效应 。 采 用 Bayesian 法 和 潜 


节 效 应 分 析 、 交 叉 滞 后 模型 中 的 纯 个 体内 的 洪 调 
节 效 应 分 析 ( 方 程 (14) 和 (15)) 、 洪 变量 增长 模型 中 
的 Sx 对 ax 一 Sy 关系 的 调节 效应 分 析 。 第 二 类 是 调 
节 变 量 Z 不 随时 间 变 化 的 调节 效应 分 析 。 如 多 层 
模型 中 的 2x(1 一 TD) 调节 效应 分 析 、 交 叉 滞 后 模型 
中 的 纯 个 体 间 的 潜 调 节 效 应 分 析 ( 方 程 (18)~ 
(20))、 洪 变量 增长 模型 中 调节 项 为 axZ Syz 的 
调节 效应 分 析 。 第 三 类 是 调节 变量 Z 随时 间 变 化 


表 2 多 层 结构 方程 模型 的 调节 分 析 的 部 分 参数 估计 结果 


参数 Yoo You Yo2 Vio 20 Ej Hoj Hy 
参数 估计 1.545 0.445 0.009 0.196 —0.052 0.167 0.057 1.002 
标准 误 0.015 0.033 0.063 0.041 0.031 0.006 0.007 0.079 


可 靠 区 间 下 限 1.513 0.384 —0.128 0.120 —0.12 0.156 0.042 0.849 
可 靠 区 间 上 限 1.571 0.513 0.119 0.276 0.003 0.178 0.070 1.154 
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的 调节 效应 分 析 。 如 多 层 模 型 中 的 1x(1 一 1) 调 节 
效应 分 析 、 交 又 滞 后 模型 中 的 三 种 潜 调 节 效 应 的 
整合 分 析 、 洪 变量 增长 模型 中 调节 项 为 axaz 或 
SxSz 的 调节 效应 分 析 。 
在 阐述 了 基于 多 层 模 型 、 多 层 结构 方程 模型 、 
交叉 滞后 模型 和 潜 变 量 增长 模型 的 多 种 纵向 数据 
的 调节 效应 分 析 后 ， 对 基于 四 类 模型 的 调节 效应 
分 析 方 法 进行 了 综合 比较 ， 接 着 总 结 出 一 个 纵向 
数据 的 调节 效应 分 析 流 程 图 ， 并 用 示例 演示 如 何 
进行 纵向 数据 的 调节 效应 分 析 。 但 是 ， 有 关 纵 向 
数据 的 调节 效应 分 析 仍 存在 不 足 ， 尚 需 进 一 步 深 
入 人 研究 和 拓展 。 
第 一 ,本 文 只 考虑 了 自 变 量 和 因 变 量 有 纵向 
数据 的 情况 。 Montoya (2019) 阐 明了 两 水 平 被 试 内 
设计 ( 仅 因 变 量 有 纵向 数据 ) 的 调节 效应 分 析 。 具 体 
地 ,研究 者 对 同一 组 被 试 施加 了 两 种 不 同 水 平 的 
实验 处 理 ( 自 变量 是 二 分 类 别 变 量 ), 调节 变量 Z 
是 不 随时 间 变 化 的 连续 变量 ， 因 变量 了 是 随时 间 
变化 的 连续 变量 (两 种 处 理 下 的 因 变 量 页 和 Yo) 
两 水 平 被 试 内 设计 的 调节 效应 分 析 的 方法 是 : 由 
Yur = 多 ~ 产生 新 变量 Yam 如 果 回 归 方 程 
Yje =dp+pbpZ+snp 的 系数 bp 显著 ,就 表示 调节 
效应 显著 。 此 外 ， 如 果 自 变量 和 调节 变量 都 是 不 
随时 间 变 化 的 连续 变量 ， 因 变量 是 随时 间 变 化 的 
连续 变量 ,可 以 使 用 多 层 模型 的 2 x (2 一 1) 调 节 效 
应 分 析 。 如 果 自 变量 是 不 随时 间 变 化 的 连续 变量 ， 
调节 变量 和 因 变 量 是 随时 间 变 化 的 连续 变量 ,可 
以 使 用 多 层 模 型 的 1 x (2 一 1) 调 节 效 应 分 析 。 
B, 本文 仅 涉及 测量 次 数 较 少 (1< 10) 的 纵 
向 数据 的 调节 效应 分 析 ， 如果 是 重复 测量 次 数 较 
多 (1 三 10) 的 纵 问 数据 ， 常 称 为 密集 追踪 数据 
(Intensive Longitudinal Data， 唐 文清 等 , 2020; 郑 舒 
方 等 ， 2021)， 如 何 进行 纵向 数据 的 调节 效应 分 析 
呢 ? 密集 追踪 数据 常用 动态 结构 方程 模型 (Dynamic 
Structural Equation Models，DSEM) 进 行 数据 分 析 
( 郑 舒 方 等 , 2021; McNeish & Hamaker, 2020)。 动 
态 结构 方程 模型 是 将 交叉 滞后 模型 和 多 层 模 型 相 
结合 ,， 其 本 质 是 在 多 层 模 型 中 添加 自 回 归 效 应 和 
滞后 效应 (方程 (31) 和 (32))。 基 于 动态 结构 方程 模 
型 的 密集 追踪 数据 的 调节 效应 分 析 可 表示 为 
层 -1: Yj = boj + By + Pej Xap + & BV) 
层 -2 : Boj = Yoo + YorZj + Ho; (32a) 


Pi; =o t MZ; + 4; (32b) 


Boj =V + 72j + Ma; (32c) 
其 中 ,系数 yio 表 示 自 回归 效应 , yoo BEAN TT BUN, 
系数 yl 表示 不 随时 间 变 化 的 层 2 调节 变量 乙 对 自 
回归 效应 Ya- Y; 的 调节 效应 ， 系 数 yo, 表示 调 
节 变 量 2 对 滞后 效应 ,1y 一 了 的 调节 效应 ， 如 
RAB ya ABO 显著 不 为 0， 就 表示 调节 效应 
显著 。 在 动态 结构 方程 模型 中 , 层 1 变量 按 组 均 
值 中 心 化 , Mplus 软件 自动 将 层 1 变量 分 解 为 个 体 
内 部 分 和 个 体 间 部 分 两 个 潜 变 量 ,2 按 总 均值 进 
行 中 心 化 (McNeish & Hamaker, 2020)。 基 于 动态 
结构 方程 模型 的 调节 效应 分 析 的 不 足 在 于 , 目前 
仅 能 分 析 不 随时 间 变 化 的 调节 变量 Z 的 调节 效应 
(Ozkok et al., in press). 

第 三 ,本 文 仅 涉 及 纵向 数据 的 调节 效应 分 析 
方法 , 已 有 研究 还 介绍 了 纵向 数据 的 中 介 效 应 分 
析 方 法 ( 刘 国 芳 等 , 2018; 方 杰 等 , 2021)， 那 么 纵 
向 数据 的 中 介 和 调节 效应 整合 在 一 起 ， 如 何 进行 
纵向 数据 的 有 调节 的 中 介 效 应 分 析 呢 7? 目前， 纵 
向 数据 的 有 调节 的 中 介 效 应 检验 可 采用 多 层 模 型 
(Bauer et al., 2006) 和 多 层 自 回归 模型 (Zhang et al., 
2018) 进 行 分 析 。 但 这 两 种 方法 都 没有 将 变量 分 为 
随时 间 变 化 的 个 体内 部 分 和 不 随时 间 变 化 的 个 体 
间 部 分 。 那 么 ,在 纵向 数据 的 有 调节 的 中 介 效 应 
分 析 中 ， 如 何 实现 随时 间 变 化 的 个 体内 效应 和 不 
随时 间 变 化 的 个 体 间 效应 的 有 效 分 离 呢 ? 还 有 待 
深入 研究 。 

纵向 数据 在 心理 学 和 其 他 社 科 领域 研究 中 经 
常 遇 到 ， 但 不 仅 纵向 数据 的 调节 分 析 的 方法 学 研 
究 沿 后 于 中 介 分 析 方 法 的 发 展 ， 而 且 纵向 数据 的 
调节 分 析 的 应 用 研究 在 国内 也 很 少 ,本 文 试 图 为 
应 用 研究 者 提供 相应 的 方法 学 支持 。 实 际 上 ，, 纵 
向 数据 的 调节 效应 的 分 析 方法 仍 处 莲 勃 发 展期 ， 
还 有 诸多 问题 冰 竺 解决。 例如， 基于 交叉 滞后 模 
型 的 纵向 数据 的 调节 效应 分 析 ， 目 前 仅 能 提供 非 
标准 化 解 ,那么 如 何 提供 标准 化 解 呢 ? 又 如 ， 纵 
向 数据 的 调节 效应 分 析 目 前 主要 通过 Mplus 软件 
实现 , 这 在 一 定 程度 上 限制 了 应 用 研究 者 进行 纵 
向 数据 的 调节 效应 分 析 ， 如 何 才能 在 其 他 软件 上 
进行 纵向 数据 的 调节 效应 分 析 呢 ?再 如 ,虽然 潜 
调节 结构 法 方程 法 优势 突出 ,但 很 多 应 用 研究 者 
放弃 该 法 的 一 个 重要 原因 就 是 该 法 耗 时 长 ， 如 何 
才能 减少 潜 调 节 结 构 方 程 法 的 耗 时 呢 ? 方法 的 进 
步 给 研究 者 提供 了 一 个 深入 理解 和 应 用 纵向 数据 
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进行 调节 效应 分 析 的 机 会 ， 相 信 随 着 纵向 数据 的 
调节 分 析 方 法 研究 的 深入 , 会 不 断 增加 我 们 对 纵 
向 数据 的 调节 效应 问题 的 理解 。 
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Moderation analysis for longitudinal data 
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of Psychology, South China Normal University, Guangzhou 510631, China) 


Abstract: At present, the analysis of moderating effect is mainly based on cross sectional data. This article 
discusses how to analyze the moderating effect with longitudinal data. If the independent variable X and the 
dependent variable Y are longitudinal data, longitudinal moderation model can be divided into three categories 
according to the type of moderator: time-invariant moderator, time-variant moderator, and moderator 
generated from X or Y. Four types of longitudinal moderation analysis approaches are summarized: 1) Multilevel 
model (MLM); 2) Multilevel structural equation model (MSEM); 3) Cross-lagged model (CLM); 4) Latent 
growth model (LGM). It is found that the decomposition of the moderating effect and the use of latent 
moderating structural equation (LMS) method are the two characteristics of the moderation analysis for 
longitudinal data. A procedure for longitudinal moderation analysis is proposed. An example is illustrated 
with the software Mplus. Directions for future research on moderation analysis of longitudinal data are 
discussed, such as the moderation analysis for intensive longitudinal data based on the dynamic structural 
equation model. 

Keywords: longitudinal data, moderation effect, multilevel model, multilevel structural equation model, cross- 
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附录 : 
基于 交叉 滞后 模型 的 纯 WxW 调节 效应 分 析 的 Mplus 程序 : 
DATA: FILE = 1.dat; 


VARIABLE: NAMES = x1-x4 ml-m4 yl-y4; 
USEVARIABLES = x1-x4 yl-y4; 
DEFINE: CENTER xl-x4 yl-y4(GRANDMEAN); ! 按 总 均值 中 心 化 
ANALYSIS: ESTIMATOR = BAYES; MEJ Bayesian 法 
TYPE = RANDOM; 


PROCESSORS = 8; ! 多 个 处 理 器 ,加 快运 算 速 度 
CHAINS = 2; ! 两 条 马尔 科 夫 链 
THIN =2; ! 抽 取 后 验 分 布 的 间隔 


FBITERATIONS = 30000;  ! 马 尔 科 夫 链 迭代 3 万 次 
MODEL: xb by xl-x4@1; NZE x 的 个 体 间 部 分 ， 即 随机 截 距 
yb by yl-y4@1; NE y 的 个 体 间 部 分 ， 即 随机 截 距 
xwl by x1@1; x1@0.01; NZE x 的 个 体内 部 分 
xw2 by x2@1; x2@0.01; 
xw3 by x3@1; x3@0.01; 
xw4 by x4@1; x4@0.01; 
ywl by yl1@1; yl@0.01; NZE y 的 个 体内 部 分 
yw2 by y2@1; y2@0.01; 
yw3 by y3@1; y3@0.01; 
yw4 by y4@1;  y4@0.01; 
xb yb with xwl-xw4@0 ywl-yw4@0; ! 随 机 截 距 和 个 体内 部 分 相关 为 0 


as 


yb on xb; ! 个 体 间 效 应 
xw4 on xw3(beta_x);  ! 自 回归 效应 


xw3 on xw2(beta_x); 


xw2 on xwl(beta_x); 

yw4 on yw3(beta_y); 

yw3 on yw2(beta_y); 

yw2 on ywl(beta_y); 

xw4 on yw3(gama_x); Hira OV 

xw3 on yw2(gama_x); 

xw2 on ywl(gama_x); 

yw4 on xw3(gama_y); 

yw3 on xw2(gama_y); 

yw2 on xwl(gama_y); 

xwl-xw4 ywl-yw4 WITH xwl-xw4@0 ywl-yw4@0; ! 限制 所 有 相关 为 0 

xwl with ywl; xw2 with yw2; ! 人 允许 自由 估计 的 相关 

xw3 with yw3; xw4 with yw4; 

xwywl | xwl XWITH ywl; ! 构 建 调节 项 
xwyw?2 | xw2 XWITH yw2; 
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xwyw3 | xw3 XWITH yw3; 
xw2 on xwywl(cx); 
xw3 on xwyw2(cx); 
xw4 on xwyw3(cx); 
yw2 on xwywl(cy); 
yw3 on xwyw2(cy); 
yw4 on xwyw3(cy); 

OUTPUT: TECH8; 

PLOT: TYPE = PLOT2; 

TE: 如 果 要 进行 跨 层次 调节 效应 分 析 ( 以 Xs;Yy1) 为 例 )， 只 需要 添加 如 下 语句 
xbywl | xb XWITH ywl; ! 构 建 调节 项 


xw2 on xbywl(clx); 


基于 多 层 结构 方程 模型 的 调节 效应 分 析 的 Mplus 程序 : 
DATA: FILE = 1.csy; 
VARIABLE: NAMES =x zy person; 
CLUSTER = person; 
ANALYSIS: TYPE = TWOLEVEL RANDOM; 
ESTIMATOR = BAYES; 
ALGORITHM = GIBBS(RW); 
POINT = MEAN;”! 用 均值 作为 点 估计 
CHAINS = 2; 
PROCESSORS = 3; 
BITERATIONS = 50000 (20000); ! 47K PAR BEIR IER 2 万 次 ， 最 多 5S 万 次 
THIN = 10; 
MODEL: %within% 
xw by x@1; x@.01; 
zw by z@1; z@.01; 
sl | y on xw; 


s2 | y on zw; 
xw with zw; 
xzw | xw XWITH zw; ! 构 建 个 体内 的 调节 项 
y on xw zw xzw; ! LMS 法 求 调节 效应 y30 
Y%between% 
xb by x@1; x@.01; 
zb by z@1; z@.01; 
sl on zb; ” ! 随 机 系数 预测 法 求 调节 效应 yi 
xzb | xb XWITH zb; 。 ! 构 建 个 体 间 的 调节 项 
y on xb zb xzb; !LMS 法 求 调节 效应 y03 
xb with zb; 
OUTPUT: CINTERVAL(hpd) TECHS; 
plot: type = plot2; 
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基于 潜 变 量 增长 模型 的 调节 效应 分 析 的 Mplus 程序 : 


DATA: 
VARIABLE: 
ANALYSIS: 


MODEL: 


OUTPUT: 


FILE = 1.dat; 
NAMES = x1-x4 z1-z4 yl-y4; 
ESTIMATOR = BAYES; 
TYPE = RANDOM; 
POINT = MEAN; 
CHAINS = 3; 
PROCESSORS = 3; 
BITERATIONS = (20000); 
BCONVERGENCE = 0.025;”! 马 尔 科 夫 链 收 敛 标准 
THIN = 10; 
IX SX | x1@0 x2@1 x3@2 x4@3; ! 构 建 LGM 模型 
IZ SZ | z1@0 z2@1 z3@2 24@3; 
IY SY | y1@0 y2@1 y3@2 y4@3; 
intl | SX XWITH SZ; ，! 构 建 变化 率 的 调节 项 
int2 | IX XWITH IZ; ! 构 建 初始 水 平 的 调节 项 
SY ON SX SZ intl; 
IY ON IX IZ int2; 
CINTERVAL(hpd) TECH8; STDYX; 


PLOT: TYPE = PLOT2; 


